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ABSTRACT 


INTRODUCTION 


The  sound  generated  by  uniform  flow  over  a  heated  cylinder 
is  calculated  for  ?  of  Reynolds  and  Grashof  numbers. 
The  time-dependent,  incompressible  flow  field  is  computed 
using  a  vorticity-stream  function  method  containing  the 
Boussinesq  approximation.  Using  the  computed  incompres¬ 
sible  flow  field,  the  radiated  sound  is  computed  using  a  low 
frequency  Green’s  function  technique.  Snapshots  of  the 
vorticity  and  temperature  fields  and  time  histories  of  the 
radiated  sound  pressure  are  presented  and  compared  with 
exi.>i.uig  results. 


NOMENCLATURE 


So 

Sound  speed 

Cd 

Drag  coefficient 

Fb 

Buoyant  force 

Gr 

Grashof  number 

Gr/Re^ 

Buoyancy  parameter 

Mo 

Mach  number 

Nu 

Nusselt  number 

P 

Pressure 

Pr 

Prandtl  number 

P. 

Sound  pressure 

r,R,  Ro,  Rmio, 

Radius 

Re 

Reynolds  number 

St 

Strouhal  number 

t 

Time 

T,ToT.^ 

Temperature 

U,V,V„Ve 

Velocity  components 

X 

Coordinate  position 

y 

Spatial  coordinate 

Y 

Potential  function 

a 

Volume  Expansion 

K 

Thermal  conductivity 

P 

Dynamic  viscosity 

D 

Vorticity 

V 

Stream  function 

P.Po 

Density 

9 

Angular  coordinate 

Transformed  coordinates 

Thit  paper  is  declared  work  of  Ihe  US.  Government  and  is 
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The  noise  radiated  from  the  low  Mach  numcci  flow 
around  a  cylinder  is  a  fundamental  acoustic  process,  and  it 
is  this  kind  of  time  dependent  wake  flow  which  can  be 
found  in  various  engineering  fields,  for  example  cooling 
fans.  The  ability  to  compute  the  sound  generated  by  a  fluid 
flow  has  been  a  long  term  goal  of  acoustics  community. 
Pioneering  work  by  Lighthill*  and  others  has  led  to  many 
approaches  for  determining  flow  generated  sound.  A  prob¬ 
lem  conuiion  to  most  approaches  is  the  assumption  that  the 
source  terms  are  either  known  from  experiment  or  can  be 
adequately  modeled. 

It  is  safe  to  say  that  the  terms  involved  in  the  com¬ 
plete  acoustic  formulation,  for  example  Lighthill’ s  acoustic 
analogy,  are  very  difficult  to  measure  and  must  be 
simplified  for  either  calculation  or  determination  by  experi¬ 
ment.  An  alternate  approach  to  specification  of  the  source 
terms  would  be  computation  of  the  flow  field  using  the  full 
Navier-Stokes  equations.  The  mature  field  of  Computa¬ 
tional  Ruid  Dynamics  (CRD)  makes  it  realistic  to  think  that 
numerical  solutions  to  the  Navier-Stokes  equations  can  be 
obtained  with  sufficient  accuracy  to  provide  the  inputs 
required  by  the  acoustic  equations.  A  sample  of  the  many 
numerical  studies  of  cylinder  flows  can  be  found  in  Refs. 
23.7.9  and  11. 

A  Computational  Acoustics  (CA)  approach  has  been 
used  by  several  researchers  to  compute  the  sound  generated 
by  low  Mach  number  flows  about  circular  cylinders.  It  is 
well  known  that  cylinders  have  been  the  object  of  countless 
experimental,  numerical  and  analytical  studies  over  the 
years.  Thus  the  phenomena  of  vortex  shedding  is  well 
understood  and  provides  a  suitable  platform  for  developing 
techniques  to  compute  the  dipole  radiation  at  low  Reynolds 
numbers.  Blake‘  provides  an  extensive  review  of  the  work 
associated  with  dipole  sound  radiated  from  cylinders  while 
the  numerical  results  of  Hardin  and  Lampkin,^  for  example, 
have  shown  good  agreement  in  the  radiated  sound  pressure 
obtained  from  numerical  solutions  and  experimental  data. 

The  above  mentioned  acoustic  computations  provide 
the  motivation  for  this  work.  Thus,  the  purpos,.  of  the 
numerical  experiments  presented  herein  is  to  predict  the 
radiated  sound  pressure  of  a  heated  cylinder  in  a  uniform 
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flow  field.  Computation  of  the  far  field  radiated  sound  is 
determined  using  a  two  step  approach.  First,  the  time  depen¬ 
dent  incompressible  flow  field  is  computed  using  CFD  tech¬ 
niques.  Second,  the  acoustic  pressure  is  computed  via  an 
integral  over  the  flow  field  using  the  low  frequency  Green’s 
function  pressure  equation  developed  by  Howe.*'®  In  the  fol¬ 
lowing  sections  the  CFD  and  CA  methods  will  be  presented 
along  with  results  of  the  flow  and  sound  computations. 

CFD  METHOD 


Equations  of  Motion 

Computational  Fluid  Dynamic  (CFD)  techniques 
will  be  used  to  determine  the  flow  field  data  required  in  the 
acousdc  calculations.  The  governing  two  dimensional 
equations  arc  continuity,  the  Navier- Stokes  equations  with  a 
buoyant  force  (using  the  Boussinesq  approximation),  the 
energy  equation  and  an  equation  of  stale.  They  are 


vV  =  o 

(1) 

Po-^=Pb-v^+pv^V 

(2) 

(3) 

p  =  pjl  -  a(T-To)] 

(4) 

Solution  of  the  above  equations,  for  the  geometry 
shown  in  Fig.  1,  makes  use  of  a  vorticity-stream  function  in 
polar  coordinates.  The  two  dimensional  velocity  com¬ 
ponents  in  terms  of  the  stream  function  V  are 


V,= 


ii^lL 

r  30 


(5) 


and  the  vorticity  is 


1  3(rV9)  1  3V, 

r  3r  r  39 


(6) 


The  equations  are  non-dimensionalized  with  respect 
to  the  cylinder  radius,  free  stream  velocity  and  a  reference 
temperature  difference.  To  resolve  flow  gradients  near  the 
cylinder  wall,  the  following  grid  stretching  transformation  is 
used 

r  =  e"^  0  -  Ttri  g=7i^e^  (7) 

After  applying  the  above  definitions,  the  flow  equations 
assume  the  form^ 


and 


3n  a(V,G)  3(VeO) 

^  3t  3^  dri 


/-  Gr  [  .  ar  3T 

Vg^smrtq— -cosrtn  — ■ 

^  2  3^  ^  3^ 

Re  ^  3ti^ 


^  3(V,T)  3(VeT) 

®  3t  3^  3n 


1 

PrRe 


3^ 


3^ 

3ti^ 


(9) 


(10) 


The  equations  are  solved  in  the  computational 
domain  and  Q<'n<Tl,...  with  the  initial  conditions 

of  V,  =  Vg  =  T  =  £T  =  0.  The  boundary  conditions  are 


3tl/ 


=  0  on  ^  =  0  and  0<ti<t|„^ 


y  =  -2  sinhTi^  sirtJlTi  and  0=0 


on  ^  =  and  {)<Ti<q^ 


For  cylinder  heating,  a  uniform  temperature  (T„^1.0)  is 
aprplied  on  the  cylinder  wall  (^=0)  after  a  i=l  startup  time. 
Periodic  boundary  conditions  are  apjplied  at  t|=1  andi)^ 
for  0<£<£_—  and  the  vorticity  SI  is  extrapx>lated  to  the 
outflow  boundary  within  the  region  0|in,  =  ±45°. 


Additionally,  several  flow  and  heat  transfer  quanti¬ 
ties  are  evaluated.  The  cylinder  surface  pressure  is  com¬ 
puted  using 

a 

P(e)=-^|n(e)t]t=<,d9  +  -^sin0  (ii) 

the  drag  coefficient  is  obtained  from 

2k 

Cd  =  f  [  n(e)t]t=o  +  n(9)6=o]  sin0  dO  ( 1 2) 


while  the  Nusselt  number  and  its  average  are  obtained  using 


Nu(0)  =  -2 


3T(9) 


5=0 


and 


(13) 


Nu;,vo=^|Nu(0)d0 


(14) 


3^'  ®  ^ 


(8) 


Solution  Method 
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The  vorticity  and  tempjerature  equations  were  solved 
using  a  method  developted  by  Torrance  and  Rockett'^  with 
modification’  to  the  body  force  terms.  Both  equations  can 


CA  METHOD 

Neglecting  the  Mach  number  dependence,  Eqn.  (17) 
Results  of  flow  field  calculations  yield  time  histories  dimensional  flow  becomes 

of  the  vorticity,  stream  function  and  velocity  components 
throughout  the  computational  domain.  These  data  will  be  = 

used  to  compute  the  sound  generated  by  the  cylinder’s  wake 

+  X2|dyi  jdyjfi  (19) 

The  method  used  herein  is  that  of  Howe^'*  who  con¬ 
sidered  the  sound  generated  by  a  Dody  in  a  non-uniform  ^  transformed  coordinates,  the  sound  integral  term  (within 

flow.  The  flow  to  be  considered  is  isentiopic  and  has  as  its  brackets)  becomes  a  function  of  (^1)).  If  the 

only  important  noise  source  a  dipole  (div  (d  x  V)).  The  considered  compact,  one  can  neglect  the  vari- 

source  is  assumed  to  be  compact,  and  the  mean  flow  has  a 

Mach  number  satisfying  M,^-«l  which  implies  that  the  flow  computed  at  each  time  step,  and  the  acoustic  pressure 

is  incompressible  with  a  constant  speed  of  sound.  Howe  determined  by  computing  the  time  derivative  of  the 

determined  that  this  kind  of  flow  can  be  approximated  by  integral.  This  was  accomplished  by  using  five  point 

the  convective  wave  equation,  and  considering  the  case  of  a  averaging  on  the  sound  integral  data  and  4th  order  central 

turbulent  eddy  convected  past  a  rigid  body.  Howe  deter-  difference  expressions  for  the  time  derivative.  Results  of  the 

mined  that,  using  his  low  frequency  Green’s  ftmetion,  the  acoustic  computations  will  be  presented  in  the  next 

sound  radiation  has  the  form  section. 


Po 

47ta„x^ 


B~-C,-^J  Hx^-V-  -vj  d’y  (16) 

where  cj  =  l/(4rcaol5f  I)  and  ?  is  the  potential  of  the 
incompressible  flow  about  the  body  and  is  of  unit  speed  at 
large  distances  from  the  body.  The  quantity  within  the 
square  brackets  is  evaluated  at  the  retarded  time 
t-  (Itl-Mo'lO/a. 

Introducing  a  velocity  which  is  parallel  to  the  nor¬ 
mal  of  the  wave  fronts  reaching  the  observation  point  5f  and 
noting  that  the  wavelength  X  of  the  sound  satisfies  the  rela¬ 
tion  X/2R  =  l/StM„,  the  conditions  for  application  of  the  low 
frequency  Green’s  function  method  have  been  met.  Thus  the 
expression  for  the  sound  radiated  to  the  far  field  observer 
located  at  t  at  time  t  is 


RESULTS 


Natural  Vortex  Shedding 


The  first  case  considered  was  the  effect  of  grid  size 
on  natural  vortex  shedding  at  three  Reynolds  numbers, 
Re=100,  200,  and  500.  As  pwinted  out  by  Rosenfeld,'^  solu¬ 
tion  convergence  may  not  occur  for  cylinder  grid  sizes  less 
than  approximately  257x257.  However,  to  keep  the  compu¬ 
tation  times  at  a  reasonable  level,  grid  densities  of  IClxlOl 
and  201x201  were  used  to  demonstrate  the  effect  of  grid 
size  on  the  solution. 


The  results  are  given  in  Figs.  2-6  for  101^  grid  and 
in  Figs.  7-13  for  the  201*  grid.  Examination  of  the  vorticity 
contours  shows  how  much  flow  detail  is  lost  using  the  101* 
grid.  At  the  same  time  comparable  results  were  obtained  for 
the  surface  vorticity  and  pressure.  Additionally,  the  com- 
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puied  drag  coefficients  are  given  in  Table  I  below  and  agree 
well  with  existing  data,  for  example  Son  and  Hanratty." 
These  data  suggest  that  reasonable  values  of  surface  data 
and  drag  can  be  obtained  using  the  smaller  grid  size. 

Table  1.  Computed  drag  coefficients  for  natural 
vortex  shedding. 


101^  Grid 

201^  Grid 

Re 

Cd 

Cd 

100 

1.095 

1.156 

200 

0.951 

1.083 

500 

0.823 

1.027 

Comparison  of  the  computed  sound  integral  shows 
substantial  changes  in  the  signatures  for  the  two  grids.  The 
101^  grid  sound  data  for  Re=200  agrees  well,  at  least  quali¬ 
tatively,  with  the  results  of  Hardin^  et  al.  The  significant 
change  in  sound  signatures  using  the  larger  grid  makes  one 
question  if  the  solution  has  converged  to  its  final  value.  It 
should  also  be  noted  that  the  estimated  Strouhal  number  for 
Re=2(X)  yields  St=0.17  as  compared  to  the  accepted  experi¬ 
mental  value  of  St=0. 19-0.20.  A  ’ower  value  of  St  is  found 
for  aU  the  data  and  reflects  what  is  believed  to  be  dissipation 
introduced  by  the  numerical  method  in  this  highly  convec¬ 
tive  flow.  However,  it  is  felt  that  the  results  for  bulk  flow 
properties  and  signature  characteristics  are  representative  of 
the  cases  considered  herein.  The  current  results  and  the  data 
of  Rosenfeld'®  indicate  the  need  to  use  the  201x201  grid  for 
the  heated  cylinder  computations  to  be  presented  in  the  fol¬ 
lowing  section. 

Uniform  Wall  Temperature 

Numerical  experiments  were  conducted  to  examine 
the  effects  of  uniform  wall  heating  and  variation  of  the 
buoyancy  parameter  Gr/Re^  on  the  vonicity  ahd  radiated 
sound.  The  results  of  these  calculations  are  presented  in 
Figs.  14-49  and  Table  2  below  for  Re=100,  200  and  500  and 
Gr/ReM>-5.0. 

After  a  startup  time  of  t=1.0,  the  unifotm  wall  tem¬ 
perature  T=1.0  was  applied.  The  isoterms  at  t=100  were 
constructed  using  a  range  of  0.5<T<1.0  to  limit  the  contours 
in  the  high  temperature  gradient  region  near  the  cylinder 
wall.  The  vorticity  and  isotherm  data  show  that  as  the  buoy¬ 
ancy  parameter  is  increased,  the  symmetrical  nature  of  the 
wake  is  lost  and  large  variations  of  vorticity  and  tempera¬ 
ture  are  found  in  the  upper  wake  region.  For  Re=100  and 
200,  the  vorticity  and  temperamre  contours  show  elimina¬ 
tion  of  the  wake  fluctuations  when  Gr/Re*=2.5  and  5.0. 
However,  the  wake  fluctuations  of  vorticity  and  temperature 
remain  for  all  values  of  Gr/Re^  for  the  Re=500  case.  In  all 
cases,  the  isotherms  follow  the  vorticity  and  heating  of  the 


fluid  only  occurs  in  the  wake  region  of  the  flow.  These 
trends  have  been  obtained  in  previous  numerical  studies^'^'^. 

The  distribution  of  surface  vorticity,  pressure  and 
Nusselt  number  have  also  been  presented  for  the  varioti'  Re 
and  Gr/Re^  at  t=100.  These  data  show  the  increase  in  mag¬ 
nitude  of  vorticity,  pressure  and  Nusselt  number  with 
increasing  Re  and  Gr/Re^.  The  data  also  show  the  largest 
variation  in  the  rear  portion  of  the  cylinder  which  is  strongly 
affected  by  the  wake  conditions.  The  computed  drag 
coefficients  and  average  Nusselt  numbers  are  given  in  Table 
2  and  show  an  increase  with  both  Re  and  Gr/Re^.  The 
Nuavo  data  agrees  well  with  the  extrapolated  data  (t— »")  of 
Jain  and  Goal^  for  Re=100  and  200. 

Table  2.  Computed  drag  coefficients  and 
Nusselt  numbers  for  various  Re  and  Gr/Re^. 


Re 

Gr/Re^ 

Cd 

Nuavo 

100 

0 

1.156 

— 

0.25 

1.170 

5.077 

0.0 

1.190 

5.115 

1.0 

1.198 

5.235 

2.5 

1.349 

5.561 

5.0 

1.622 

6.006 

200 

0 

1.083 

0.25 

1.095 

7.355 

0.5 

1.098 

7.375 

l.O 

1.160 

7.590 

2.5 

1.179 

7.790 

5.0 

1.549 

8.313 

500 

0 

1.027 

0.25 

1.074 

11.867 

0.5 

1.034 

12.228 

1.0 

1.032 

12.252 

2.5 

1.188 

12.811 

5.0 

1.589 

12.642 

Using  Eqn.  (19),  the  time  history  of  the  radiated 
sound  integral  (proportional  to  the  radiated  sound  pressure) 
was  computed  for  several  angular  positions,  (O°<0<9O“),  for 
each  Re  and  Gr/Re^  Examination  of  the  data  shows  the 
periodic  nature  of  the  acoustic  radiation  and  the  increase  in 
magnitude  with  increasing  Reynolds  number  and  the  change 
in  signature  character  for  increasing  Gr/Re^  Although  not 
shown,  the  radiated  sound  becomes  =  zero  for 
Gr/Re^2.5  and  5.0  for  Re=100  and  Gr/Re^5.0  for  Re=200 
which  reflects  the  elimination  of  vorticity  fluctuations  in  the 
heated  wake.  Similar  results  were  obtained  by  Mautner’  for 
cylinder  wake  flows  in  stratified  flow.  TTie  details  of  the  sig¬ 
natures  can  be  obtained  from  the  figures.  However,  it  should 
be  noted  that  as  the  wall  heating  and  buoyancy  begin  to 
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modify  tl»c  wi<ke  vorticity,  the  radiated  sound  signatures 
loose  the  sirtusoidal  nanire  found  for  T=0  and  reflect  the 
wake  fluctuations  over  O<0<9O°.  This  would  suggest  that 
not  only  is  the  magnitude  of  the  radiated  sound  increased  in 
the  heated  wake,  but  a  larger  spectrum  of  frequencies  of 
sound,  in  addition  to  the  fundamental,  will  be  radiated  into 
the  far  field. 


CONCLUSIONS 


The  problem  of  laminar  mixed  convection  firom  a 
heated  cylinder  has  been  investigated  numerically  by  solv¬ 
ing  the  Navier-Stokes  and  energy  equations  over  a  range  of 
Re  and  Gr/Re^.  The  radiated  sound  has  been  computed 
using  a  low  frequency  Green’s  function  method.  The  effect 
of  grid  size  on  solution  convergence  dictated  the  use  of  a 
201^  grid  for  the  heated  cylinder  flow  and  acoustic  computa¬ 
tions.  The  natural  shedding  results  were  in  good  agreement 
with  previous  results  as  were  the  trends  in  the  temperature 
contours  for  various  values  of  the  buoyancy  parameter.  The 
radiated  sound  results  were  in  qualitative  agreement  with 
existing  data  for  Re=200  using  the  101^  grid.  However,  low 
values  of  Strouhal  number  indicate  excessive  dissipation 
introduced  by  the  numerical  method.  The  acoustic  results 
show  significant  changes  in  the  signature  shape  with 
increasing  Gr/Re^  and  magnitude  increases  with  increasing 
Re.  The  computed  values  of  drag  coefficient,  surface  pres¬ 
sure,  vorticity  and  Nusselt  number  agree  well  with  existing 
data.  Even  though  the  small  grid  size,  (101^),  allows  for  rea¬ 
sonable  values  of  bulk  parameters,  for  example  drag,  higher 
density  grids  are  required  to  obtain  solution  convergence. 
The  change  in  character  of  the  radiated  sound  from  the  101^ 
grid  to  the  201^  grid  suggests  that  a  finer  grid  may  be 
required  to  obtain  solution  convergence.  Finally,  the  natuie 
of  the  radiated  sound  in  the  heated  wake  suggests  that  a 
much  broader  range  of  frequencies  will  be  present  in  the  far 
field  sound. 
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Fig.  1 .  Geometry 
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Fig  5.  Time  history  of  the  sound  integral 
for  the  101x101  grid  and  Re=100. 

Fig-  2.  Vorticity  contours  for  101x101  grid. 
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Fig.  3.  Cylinder  surface  voriicity  distributions 
for  the  101x101  grid. 
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Fig.  4.  Cylinder  surface  pre.ssurc  distributions 
for  the  lOlx  101  grid. 


Fig  6.  Time  history  of  the  sound  integral 
for  the  101x101  grid  and  Re=200. 
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Fig  7.  Time  history  of  the  sound  integral 
for  the  101x101  grid  and  Rc=500 
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Fig.  8.  Vorticity  contours  for  201x201  grid. 
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Fig  1 1 .  Time  history  of  the  sound  integral 
for  the  201x201  grid  and  Re=100. 
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Fig.  9.  Cylinder  surface  vorticity  distributions 
for  the  201x201  grid. 


Fig  12,  Time  history  of  the  sound  integral 
for  the  201x201  grid  and  Re=200. 
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Fig.  10.  Cylinder  surface  pressure  distributions 
for  the  201x201  grid. 
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Fig  13.  Time  history  of  liic  sound  inicgr.il 
for  llic  201  .n20l  grid  andRc-.SOO 
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Fig.  14.  Vorticity  contours  and  isothcmrs 
for  Re=100andGr/Re^0.25. 
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Fig.  15.  Vorticity  contourr  and  isotherms 
for  Re=100  and  Gr/Re^=0.5. 


Fig.  16.  Vorticity  contours  and  isotherms 
for  Re=100  and  Gr/Re^l  .O. 


Fig  17.  Vorticity  contours  and  isotherms 
lor  Rc=lOO  and  Gr/Re^2.5. 


Fig.  18.  Vorticity  contours  and  isotherms 
for  Re^^lOO  and  Gr/Re^5.0. 


Fig.  19,  Cylinder  surface  vorticity 
for  Rc=  100  and  Gr/Rc^.=0-5.0. 
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Fig  20  Cylinder  surface  prr^ssuic 
for  Rej-  100  and  Gr/Rc^O  5  0 
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Fig.  21.  Cylinder  surface  NusscU  number 
for  Rej=100  and  Gr/Re^O.25-5.0. 
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Fig.  23.  Time  history  of  the  sound  integral 
for  Re=100  and  Gr/Re^^.5. 
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Fig.  25.  Vorticily  contours  and  isotherms 
for  Rc"200  and  Gr/Rc^O.25. 


Fig.  26.  Vonicity  contours  and  isoilicnn.s 
for  Rc=200  and  Cr/Rc^  ()..5. 


rig.  28.  Vonicity  contours  i 
for  Re=200  and  Gr/Re^2.5. 
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Fig.  34.  Time  history  of  the  sound  integral  Fig.  35.  Time  history  of  the  sound  integral  Fig.  36.  Time  history  of  the  sound  integral 

for  Re=200  and  Gr/Re*=0.5.  for  Re=200  and  Gr/Re^l.O.  for  Re=200  and  Gr/Re^2  5 
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Fig.  37.  Vorticity  contours  and  isotherms 
for  Re=500  and  Gr/Re^=0.25. 
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Fig.  38.  Vorticity  contours  and  isotherms 
for  Re=500  and  Gr/Re^=0.5. 
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Fig.  39.  Vorticity  contours  and  isotherms 
for  Rc=500  and  Gr/Re^I  .O. 
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a)  Vorlidty 


Fig.  40.  Votticity  contours  and  isotherms 
for  Re=500  and  Gr/Re’i2.5. 


b)  Temperature 

Fig.  41.  Vorticity  contours  and  isotherms 
for  Re=500  and  Gr/Re^5.0. 
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Fig.  42.  Cylinder  surface  vorticity 
for  Re=500  and  Gr/ReM)-5.0. 
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Fig.  43.  Cylinder  surface  pressure 
for  Re=500  and  Gr/ReM)-5.0. 
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Fig.  44.  Cylinder  surface  Nussclt  number 
for  Re=500  and  Gr/Re^.25-5.0. 
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Fig.  45.  Time  history  of  the  sound  integral 
for  Re=500  and  Gr/Re^=0.25. 
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Fig.  46.  Time  history  of  the  sound  integral 
for  Re=500  and  Gr/Re^=0J. 
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Fig.  47.  Time  history  of  the  sound  integral 
for  Re=500  and  Gr/Re^l  .0. 


Fig.  48.  Time  history  of  the  sound  integral 
for  Re=500  and  Gr/Re^2,5. 


Fig.  49.  Time  history  of  the  sound  integral 
for  Re=500  and  Gr/Re^5.0. 
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